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Abstract 

The shear and the bulk viscosities of the hadron gas at low temperatures are studied in the 
model with constant elastic cross sections being relativistic generalization of the hard spheres 
model. One effective radius r = 0.4 fm is chosen for all elastic collisions. Only elastic collisions 
are considered which are supposed to be dominant at temperatures T < 120 — 140 MeV. The 
calculations are done in the framework of the Boltzmann equation with the Boltzmann statistics 
distribution functions and the ideal gas equation of state. The applicability of these approximations 
is discussed. It's found that the bulk viscosity of the hadron gas is much larger than the bulk 
viscosity of the pion gas while the shear viscosity is found to be less sensitive to the mass spectrum 
' (— I '■ of hadrons. The constant cross sections and the Boltzmann statistics approximation allows one not 

Q-1 only to conduct precise numerical calculations of transport coefficients in the hadron gas but also 

to obtain some relatively simple relativistic analytical closed-form expressions. Namely, the correct 
single-component first-order shear viscosity coefficient is found. The single-component first-order 
nonequilibrium distribution function, some analytical results for the binary mixture and expressions 
for mean collision rates, mean free paths and times are presented. Comparison with some previous 
calculations for the hadron gas and the pion gas is done too. This paper is the first step towards 
calculations with inelastic processes included. 
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I. INTRODUCTION 



The bulk and the shear viscosity coefficients are transport coefficients which enter in 
hydro dynamical equations and thus are important for studying of nonequilibrium evolution 
of any thermodynamic system. 

Although experimental value of the shear viscosity was found to be negligibly small in 
the hadron gas [l|, ^ there are two more additional reasons to study the shear viscosity. 
The ffist one is the experimentally observed minimum of the ratio of the shear viscosity to 
the entropy density rj/s near the liquid-gas phase transition for different substances which 
may help in studying of the quantum chromodynamics phase diagram and finding of the 
location of the critical point [3|, [2|. For a counterexample see {4J and references therein. 
The second reason is calculation of the 7]/s in strongly interacting systems, preferably real 
ones, for testing of the conjectured lowest bound r]/s = jsj, found in some conformal field 
theories having dual gravity theories, and search of a new one. The bound 'n/s = ^ was 
violated with different counterexamples. For some reasonable ones see Also see the 

recent review j8|. The bulk viscosity may be not negligibly small and is expected to have a 



maximum near the phase transition [9|, [10|, [11 



Whether one uses the Kubo formula or the Boltzmann equation one faces nearly the same 
integral equation for transport coefficients [l2|, [l3[, [13] • The preferable way to solve it is 
the variational (or Ritz) method. Due to its complexity often the relaxation time approxi- 
mation is used in the framework of the Boltzmann equation. Though this approximation is 
inaccurate, does not allow to control precision of approximation and can potentially lead to 
large deviations. The main difficulty in the variational method is in calculation of collision 
integrals. To calculate any transport coefficient in the lowest approximation in a mixture 
like the hadron gas with very large number of components A^' one would need to calculate 
roughly A^'^ special 12-dimensional integrals if only elastic collisions are considered. Fortu- 
nately it's possible to simplify these integrals considerably and perform these calculations 
in reasonable time. 

This paper contains calculations of the shear and the bulk viscosity coefficients for the 
hadron-resonance gas in the constant elastic differential cross sections model. The calcula- 
tions are done in the framework of the Boltzmann equation with classical Boltzmann statis- 
tics, without medium effects and with the ideal gas equation of state. The constant cross 
sections and the Boltzmann statistics approximation allows one to obtain some relatively 
simple analytical closed-form expressions. 

The numerical calculations for the hadron-resonance gas presented in this paper are the 
ffist step towards calculations with number- changing (inelastic) processes taken into account. 
Then, possibly, more realistic model of strong interactions and quantum statistics will be 
implemented. These calculations can be considered as quite precise at low temperatures 
where elastic collisions dominate and equation of state is close to the ideal gas equation of 
state. 

For comparison calculations for the pion gas are performed and the results are found to 
be close to the results in [l5| where only elastic processes were considered and the ideal gas 
equation of state was used or to the results in [l6j (at zero chemical potential) where the 
Bose statistics was used instead of the Boltzmann one. The comparison is shown in fig. [H 
The discrepancies arise because of the constant cross sections approximation used in this 
paper. Smaller values of the viscosities in [l5[ for some temperature range originate because 
of the p-resonance contribution to the quasielastic vrvr cross section and the maximum of the 
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bulk viscosity is expected to be shifted towards lower temperatures. In these calculations 
the nonvanishing value of the bulk viscosity is obtained only due to the mass of pions.^ 
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FIG. 1: Dashed lines designate calculations for the pion gas performed in the constant cross sections 
approximation with the hard core radius r = 0.4 fm. The solid black line designates the 3-rd order 
bulk viscosity and the solid red line designates the 1-st order shear viscosity of the pion gas from 



15l | , where the quasielastic (with contribution of the p-resonance) semi-empirical nir cross section 
was used. 



In several papers the bulk viscosity ^ was calculated for the pion gas using the chi- 
ral perturbation theory and some other approaches with quite large discrepancies between 
quantitative results. In [19] calculations were done by the Kubo formula in rough approxima- 
tion. There the number-changing 2 ■v^ A processes were neglected too and the non- vanishing 
value of the bulk viscosity is obtained due to trace anomaly and the mass of pions. At small 
temperatures where effects of trace anomaly are small the magnitude of the bulk viscosity is 
large in compare to the results of this paper and [l^, jl6|. For example, at T = 60 MeV the 
bulk viscosity is roughly 20 times more. The calculations in [20[ are done in the framework 
of the Boltzmann equation and have divergent dependence of the ^ (T) for T — )■ because of 
weak 2 4-7-4 number- changing processes taken into account^ and at T = 140 MeV the bulk 
viscosity is nearly 44 times larger than the bulk viscosity calculated in this paper. Though 
calculations in ^0| were done in the first order approximation which does not even take 
into account elastic collisions so that this discrepancy can become smaller few times after 



^ The bulk viscosity vanishes identically in scale invariant theories 17|. Generally speaking it vanishes 
identically if the distribution functions are scale invariant [isj which explains the case of the massive 
monoatomic gas in the nonrelativistic theory [18] (where e — 3P 7^ 0) with the scale covariant energy 
spectrum 

^ This dependence should change at low temperatures, where number-changing processes are expected to 
be suppressed. See Sec. |llI]for more details. 
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more accurate calculations. The function ^ (T) may turn out to be not smooth at the middle 



temperatures and the smooth function ^ (T) is to be obtained through extrapolation. In [21 



the bulk viscosity was calculated in the framework of the Boltzmann equation with the ideal 
gas equation of state and only elastic collisions taken into account. The Inverse Amplitude 
Method was used to get the scattering amplitudes of pions. The quantitative results are 
close to the results in [l5|, ll6|. Calculations in the paper (22| are done in the framework of 
the Boltzmann equation for massless pions. There the bulk viscosity increases rapidly so 
that the ratio ^/s increases with the temperature. 

Calculation of the shear viscosity in the hadron gas, done in 23| using some approximate 
formula and in using relaxation time approximation, are in good agreement with cal- 
culations in this paper. The bulk viscosity of the hadron gas, calculated in [23] has similar 
dependence on the temperature though it is roughly 10-30 times more than the bulk vis- 
cosity calculated in this paper. At low temperature T = 100 MeV and vanishing chemical 



potentials it is 20 times more. In [25| calculation of the bulk viscosity was done for the 
hadron-resonance gas with masses less than 2 GeV using some special formula obtained 
though an ansatz for the small frequency limit of the spectral density at zero spatial mo- 
mentum [9i]. There the bulk viscosity to entropy density ratio ^/s has similar behavior to 
the ^/s in this paper and is roughly 2 times more. 

The bulk and the shear viscosities were calculated in the linear u-model with number- 



changing processes taken into account and nonideal gas equation of state [11|. This model 
has chiral symmetry restoration phase transition and can qualitatively describe pion gas. 
These calculations demonstrate example that the ratio of the shear viscosity to the entropy 
density r]/s has the minimum near the phase transition and the ratio of the bulk viscosity 
to the entropy density ^/s can have a maximum near the phase transition for some values 
of the vacuum sigma mass if the peak in the ^ is sharp enough. 

The structure of the paper is the following. Sec. [TTl contains the constant cross sections 
model description. The applicability of the used approximations is discussed in Sec. Illli 
The Boltzmann equation, its solution and general form of the transport coefficients can be 
found in Sec. IIVI The numerical calculations for the hadron gas are presented in Sec. |Vl 
In Sec. IVI Al analytical results for the single-component gas are presented. In particular 
analytical exp ression for the first order single- component shear viscosity coefficient, found 
before in |26|, is corrected while the bulk viscosity coefficient remains the same. Also the 
nonequilibrium distribution function in the same approximation is written. Some analytical 
results for the binary mixture are considered in Sec. IVI B[ Integrals of source terms needed 
for calculation of the transport coefficients can be found in Appendix |Al The general entropy 
density formula used in numerical calculations for the hadron gas can be found in Appendix 
[BI Transformations of the collision brackets being the 12-dimensional integrals which enter in 
transport coefficients can be found in Appendix O The closed-form expressions for collision 
rates, mean free paths and mean free times are included in Appendix [Dl 



II. THE HARD CORE INTERACTION MODEL 

In non-relativistic classical theory of particle interactions there is a widespread model 
called the hard core repulsion model or model of hard spheres with some radius r. For 
applications for the high-energy nuclear collisions see and references therein. The differ- 



ential scattering cross section for this model can be inferred from the problem of scattering 



of point particle on the spherical potential U{r) = oo if r < a and U{r) = if r > a [27 . 
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In this model the differential cross section is equal to a^/4. To apply this result to the 
gas of hard spheres with radius r one can notice that the scattering of two spheres can be 
considered as the scattering of the point particle on the sphere of the radius 2r, so that 
one should take a = 2r. The full cross section is obtained after integration over angles of 
the r'^dQ which results in the Anr"^. For collisions of hard spheres of different radiuses one 
should take a = ri + r2 or replace r on ^^^''^ . The relativistic generalization of this model is 
the constant (not dependent on the scattering energy and angle) differential cross sections 
model. 

The hard spheres model is classical and connection of its cross sections to cross sections 
calculated in any quantum theory is needed. For particles having spin the differential cross 
sections averaged over the initial spin states and summed over the final ones will be used.^ 
If colliding particles are identical and their differential cross section is integrated over the 
momentums (or the spatial angle to get the total cross section) then it should be multiplied 
on the factor ^ to cancel double counting of the momentum states. These factors are 
exactly the factors 'Jm next to the collision integrals in the Boltzmann equations f l4.19p . The 
differential cross sections times these factors will be called the classical differential cross 
sections. 



From the analysis of semi-empirical total elastic cross sections in [15| for tttt, ttK, ttN, 
KN and A^A^ collisions one can find that the elastic part without resonance contribution of 
the total cross sections lies approximately in the range 10 — 30 mh = 1 — 3 fm? (except for 
A^A^ cross section reaching 50 mh for small energies). From comparison of these values with 
the total cross sections in the hard spheres model one finds for the radius r = 0.28 — 0.49 fm. 
For simplification the mean value is chosen equal to r = 0.4 fm in numerical calculations 
for all hadron and resonance elastic cross sections. Then classical differential cross sections 
become equal to some one constant value too. In this approximation the cross section will 
enter in the transport coefficients as one factor. 



III. CONDITIONS OF APPLICABILITY 

First the applicability of the Boltzmann equation and of calculations of transport coeffi- 
cients should be clarified. 

Although the Boltzmann equation is valid for any perturbations of the distribution func- 
tion it should be a slowly varying function of space-time coordinate to justify that it can be 
considered as a function of macroscopic quantities like temperature, chemical potential or 
flow velocity or in other words that one can apply thermodynamics locally. Then one can 
make expansion over independent gradients of thermodynamic functions and flow velocity 
(the Chapman- Enskog method), which vanish in equilibrium. Smallness of these perturba- 
tions of the distribution functions in compare to theirs main parts ensures the validity of 
this expansion and that the gradients are small. ^ Because these perturbations are inversely 



^ It's assumed that particle numbers of the same species but with different spin states are equal. If this were 
not so then in approximation, in which spin interactions are neglected and probabilities to have certain 
spin states are equal, the numbers of the particles with different spin states would be approximately equal 
in the mean free time. With equal particle numbers theirs distribution functions are equal too. This 
allows one to use summed over the final states cross sections in the Boltzmann equation. 

^ The magnitudes of thermodynamic quantities can also be restricted by this condition or, conversely, not 
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proportional to coupling constants one can say that they are proportional to some product 
of particles' mean free paths and the gradients. So that in other words the mean free paths 
should be much smaller than the characteristic lengths on which macroscopic quantities 
change considerably^. 

The system with number- changing processes should be treated additionally. Number- 
changing processes are very important for the bulk viscosity. If the coupling constant of 
inelastic number- changing processes is small then the bulk viscosity is dominated by the 



inelastic processes [12|. Formally tending coupling constants of inelastic processes to zero the 
bulk viscosity diverges together with the nonequilibrium perturbations of the distribution 
functions, which should be small. Though this dependence can become invalid earlier if 
inelastic processes are negligible in some sense, because the limiting case with completely 
switched off inelastic processes can be well defined. That's why there is a need to specify 
reasonable conditions when inelastic processes can be neglected. Some simple ones of them 
are proposed below. It is smallness of the time t on which the system's evolution is considered 
in compare to the mean free time for number-changing processes of the k'-th^ particle species, 

tp'^ dnis]), 

^ < 1, (3.1) 



j-inel 



where is the rate of inelastic processes per particle of the k'-th species. Similar condi- 
tion, stating that the number of reactions is smaller than unity, can be used: 

tv ^S<1' (3-2) 

ng all channels 

where V is the system's volume and ^ii channels -^fc"n number of the inelastic reac- 

tions of particles of the k'-th species over all channels per unit time per unit volume. Though 
in mixtures with mean free times of inelastic processes close to each other one might need to 
use relevant sums of reaction rates instead of reaction rates for certain species. Also one can 
consider natural time-scales like relaxation times of the gradients or thermodynamic func- 
tions. Estimations of the thermal and chemical relaxation times for pions were done in j28| . 
Basing on these results one can expect that the approximate temperature where inelastic 
processes cease to be efficient is T = 120 — 140 MeV for vanishing chemical potentials. For 
nonzero chemical potentials this temperature is expected to be smaller. 

Errors due to the Boltzmann statistics used instead of the Bose or the Fermi ones were 
found to be small for vanishing chemical potentials.'' According to calculations for the 
pion gas in fl6| the bulk viscosity becomes 25% larger at T = 120 MeV and 33% larger 
at T = 200 MeV for vanishing chemical potential. Although the relative deviations of the 



restricted even if transport coefficients diverge. See also Sec. IVI Al of this paper. The smallness of the 
shear and the bulk viscosity gradients can also be checked by the condition of smallness of the T^^^'^'' 
in compare to the T^")'''" (gH). 
^ It's clear that the mean free paths should be smaller than the system's size too. 

^ Primed indexes run over particle species without regard to their spin states. This assignment is clarified 
more below. 

^ It should be mentioned that if particles of the fc-th particle species are bosons and if iik{x^) > mk then 
there is (local) Bose-Einstein condensation for them which should be treated in a special way. 
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thermodynamic quantities of the pion gas at nonvanishing chemical potential /i = 100 MeV 
are not more than 20%^ the bulk viscosity becomes up to 2.5 times more. The shear viscosity 
becomes 15% less at T = 120 MeV and 25% less at T = 200 MeV for vanishing chemical po- 
tential and 33% less at T = 120 MeV and 67% less at T = 200 MeV for the /z = 100 MeV. 
The corrections to the bulk viscosity of the fermion gas, according to calculations of the bulk 
viscosity source term not presented in this paper, are of the opposite sign and approximately 
of the same magnitude. 

The condition of applicability of the ideal gas equation of state is controlled by the di- 
mensionless parameter vn which appears in the first correction from the binary collisions in 
the virial expansion and should be small. Here v = 167rr'^/3 is the so called excluded volume 
parameter and 1/n is the mean volume per particle. One finds vn ^ 0.09 at T = 120 MeV 
and vn ^ 0.2 at T = 140 MeV for vanishing chemical potentials. Though even small correc- 
tions to thermodynamic quantities can result in large corrections for the bulk viscosity as it 
was found for the quantum statistics corrections. Along the freeze-out line (its parameters 
can be found in j23|) the vn grows from 0.07 to 0.49 with the temperature. 

One more important requirement which one needs to justify the Boltzmann equation 
approach is that the mean free time should be much larger than h/fl {fl is the character- 
istic single-particle energy) [29| or de Broglie wavelength should be much smaller than the 
mean free path [13] to distinguish independent acts of collisions and for particles to have 
well-defined on-shell energy and momentum. This condition gets badly satisfied for high 
temperatures or densities. The mean free path of particle species k' is given by the formula 
(]D16p or the formula (IDlip if inelastic processes can be neglected. The wavelength can be 
written as A^' = -^p^; where the averaged modulus of momentum of k'-th species {\pk'\) is 

/r |\ - I d'Pk'\Pk'\fi>\pk') _ 2e-'^'T{3 + 3zk> + zl,) _ j ^myT K^,2{zk>) 

where zy = nik'/T, K2{x) is the modified Bessel function of the second kind. As it follows 
from the (13. 3p the largest wavelength is for the lightest particles, vr-mesons. The elastic 
collision mean free paths are close to each other for all particle species. Hence, the smallest 
value of the ratio Xk'/^k' is for vr-mesons. Its value is close to the value of the vn and is 
exponentially suppressed for small temperatures too. At temperature T = 140 MeV and 
vanishing chemical potentials this ratio is equal to 0.18. Along the freeze-out line it grows 
from 0.12 to 0.37 with the temperature. 

If the last requirement is not satisfied one can use the Kubo formulas, for instance. In 
13 it was shown that intermediate integral equations for the shear and the bulk viscosities 



coming after linearization of the Boltzmann equation and expansion over gradients of the 
flow velocity can be reproduced from the Kubo formulas in scalar g(j)^ + theory. The 
difference is in that one has to replace particles' masses with temperature dependent ones 
and to use thermal matrix elements for elastic and inelastic collisions. Basing on this result 
an example of effective kinetic theory of quasiparticle excitations valid for all temperatures 



was presented in [13|. For further developments see [IJ], |30|, [31[. For other approaches see 



3^ . (33| . (33 | and 35] with references therein 



The relative deviations of the thermodynamic quantities grow with the temperature for some fixed value 
of the chemical potential and tend to some constant. 



7 



IV. DETAILS OF CALCULATIONS 



The Boltzmann equation and its solution 



Calculations in this paper go close to the ones in [36[ though with some differences and 
generalizations. Let's start from some definitions. Multi-indices k,l,m,n will be used to 
denote particle species with certain spin states. Indexes k', I', m', n' will be used to denote 
particle species without regard to their spin states (and run from 1 to the number of particle 
species N') and a, 6 to denote conserved quantum numbers^. Because nothing depends on 
spin variables one has for every sum over the multi-indexes 

Y.... = Y.gk'.... (4.1) 

k k' 

where Qk' is the spin degeneracy factor. The following assignments will be used 

E\ ^ \ ^ n-k Uk' ria 

rik = } nk', na = }Qaknk, Xk = — , Xk' = — , Xa = — , 
^ ^ ^ fin tl fl 



n 

k k' k 



^J'k — rp ^ IJ'a — rp ^ ^k — rp } k — rp ^ 'k — ' l^'-^J 

where Qak denotes values of conserved quantum numbers of the a-th kind of the fc-th particle 
species. Everywhere the particle number densities are summed the spin degeneracy factor 
gy appears and then gets absorbed into the ny or the by definition. All other quantities 
with primed and unprimed indexes don't differ, except for rates, mean free times and mean 
free paths defined in Appendix [D], the 7^; commented below, the coefficients A^f^/, C^f;; and, 
of course, quantities whose free indexes set indexes of the particle number densities n^. Also 
the assignment j = j^^ will be used somewhere for compactness. 
The particle number flows are^" 



T--Y.I^PiPU. (4.4) 



where the assignment fk{Pk) = fk is introduced. The energy- momentum tensor is 

d^Pk 

The local equilibrium distribution functions are 

where fik is the chemical potential of the k-th particle species, T is the temperature and 
Ufj_ is the relativistic flow 4-velocity such that U^U^ = 1 (with frequently used consequence 



In systems with only elastic collisions each particle species have their own " conserved quantum number" , 
equal to 1. 

The +,—,—,— metric signature is used throughout the paper. 
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U^dyU^ = 0). The local equilibrium is considered as perturbations of independent thermo- 
dynamic variables and the flow velocity over a global equilibrium such that they can depend 
on the space-time coordinate x'^. Though such perturbations are not the general ones and 
do not take into account all possible deviations from a chemical equilibrium. For numeri- 
cal calculations along the freeze-out line, where such deviations are important, saturation 



7-factors are used, see [37j, [38| and references therein. The chemical equilibrium implies 
that the particle number densities are equal to their global equilibrium values. The global 
equilibrium is called the time-independent stationary state with the maximal entropy. The 
global equilibrium of isolated system can be found by variation of the total nonequilibrium 
entropy functional [soj over the distribution function with condition of the total energy and 
the total net charges conservation: 

^1/1 - E / - - A) - E / - E^«^-/ ^/^. («) 

where /3, Xa are Lagrange coefficients. Equating the first variation to zero one easily gets the 
function (ES]) with W = (1, 0, 0, 0), (3 = ^ and 



^Qakf^a, (4.7) 



where fia = Xa are the independent chemical potentials coupled to conserved net charges. 

With fk = substituted in the (HJ!) and the g3!) one gets the leading contribution in 
the gradients expansion of the particle number fiow and the energy-momentum tensor 

where the projector 

A'"' = g^''' - u^U", (4.10) 
is introduced. The is the ideal gas particle number density 

nu = U,Nf^^ = ^ThlK^{z,)e^\ (4.11) 

the e is the ideal gas energy density 

, . r;^f;j,«... . ^ j ^pO;..) = ^ ^ ™.|||| - T. (4.12) 

and the P is the ideal gas pressure 
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Also the following assignments are used 

e = - = XkCk, hk = ek + T, h = — — = Xkhk, 
n ^ — ^ n ^ — ^ 

k k 

Ck = — = Zk -^ , . - 1, e = -, hk = — = Zk , . , h = -. (4.14) 

Above h is the enthalpy per particle, e is the energy per particle and hk, are the enthalpy 
and the energy per particle of the /c-th particle species correspondingly which are well defined 
in the ideal gas. 

In relativistic hydrodynamics the flow velocity needs some extended definition in 
relation to the thermodynamic quantities. The most convenient condition applied to 
is the Landau-Lifshitz condition j4o|. This condition states that in the local rest frame 
(where the flow velocity is zero though its gradient can have nonzero value) each fluid cell 
should have zero momentum and its energy and net charge densities should be related to 
other thermodynamic quantities through the equilibrium thermodynamic relations (without 
contribution of nonequilibrium dissipations). Its covariant mathematical formulation is 

(T^^ - T(°)'^'')[/^ = 0, {N^ - Nl^^>')U^ = Q. (4.15) 

The next to leading correction over the gradients expansion to the T^'^ can be written as 
expansion over the 1-st order Lorentz covariant gradients which are rotationally and space 
inversion invariant and satisfy the Landau-Lifshitz condition^^ (14.151) : 

^ 2r]WlF + ^A'^'VpUP = r] (^A^A^ + A^A^ - ^ A^'^Ap^^ Vf/^ + ^A^'^VpU'', 

(4.16) 

where for any tensor a^^^y the symmetrized traceless tensor assignment is introduced 

The equation (I4.16P is the definition of the shear rj and the bulk ^ viscosity coefficients. 
The ^A^^'^'V pW term in the fl4.16p can be considered as nonequilibrium contribution to the 
pressure which enters in the (14. 9p . 

By means of the projector fl4.10p one can split the space-time derivative (9^ as 

d,, = U^Wd, + Aid, = UpD + \/p. (4.18) 

The D becomes the time derivative in the local rest frame (where = (1, 0, 0, 0)) or in 
the nonrelativistic limit and the becomes the spacial derivative. Then the Boltzmann 
equations can be written in the form 

P'M = {p'kUpD + p^,Vp)h = C^kVk] + CrVk], (4.19) 



Also this form of r(i)^'^ respects the second law of therraodynamics [40]. Implementation of the Eckart 
condition would result in different form of the T'^^^^^^ [Sal- 
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where C^"'^'[/fc] represents the inelastic or number-changing colhsion integrals (it is dropped 
in calculations in this paper if the opposite is not stated explicitly) and Cf [fk] is the elastic 
2 -H- 2 collision integral. The collision integral Cf[fk\ has the form of the sum of positive 
gain terms and negative loss terms. Its explicit form is (cf. [l2[ 14 1) 



ct[h 



(Ppii (iVfc d^Pii 



Wkf'll ~ fkfll) 



\MMn27rmp',+p[i-pk-pu) 



(4.20) 



where jki = | if ^ ^ denote the same particle species without regard to spin states and 
•jki = 1 otherwise, \M.ki{p'ky p'w Pk, Pii)\'^ = |-^fc«P is the square of dimensionless elastic scat- 
tering amplitude, averaged over the initial spin states and summed over the final ones. Index 
1 designates that pk and pifc are different variables. Introducing Wm = Wki{p'k, Pu', Pk, Pii) as 



W, 



kl 



647r2 



^ {P'k+P'll-Pk-Pll), 



(4.21) 



one can rewrite the collision integral fl4.20p in the form as in [36 



Ct[fk] = (2vr) 



'J2lki / 

I Jpii,p'k,p'i 



fk fll 



fk fl 



(27r)3 (27r)3 (27r)3 (27r) 



W, 



kl- 



The Wki is related to the elastic differential cross section aki as [3 

Wki = s(Tki5^{p'k + p\i -pk- Pii), 



(4.22) 



(4.23) 



where s = {pk + pu)"^ is the usual Mandelstam variable. The Wki has properties 
Wki{p'k,Pu;Pk,Pii) = Wki{pk,Pii;p'k,Pii) = Wik{p'ii,p'k;Pii,Pk) (due to time reversibility and 
a freedom of relabelling of order numbers of particles taking part in reaction). And e.g. 
^kiip'k: Pw Pk, Pii) 7^ Wki{p'ii, p'k', Pii, Pk) in general case. Elastic collision integrals have im- 
portant properties which one can easily prove (36| : 



d^Pk 
(2vr)V, 



cnfk] = o, 



(4.24) 



E 



^^-QPkCklfk 



(2^)V, 



0. 



(4.25) 



Also the C^'[/fe] vanishes if fk = fjf^. 

The distribution functions fk solving the system of the Boltzmann equations approxi- 
mately are sought in the form 



p(i) 



(0) 



ciO), 



(4.26) 



The factor jki cancels double counting in integration over momentums of identical particles. The factor 
^ comes from the relativistic normalization of scattering amplitudes. 
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where it's assumed that fk depend on the x'^ entirely through the T, jj^k, or their space- 
time derivatives. Also it is assumed that \(pk\ ^ 1- After substitution of fk = f'^'' in the 
fl4.19p the r.h.s. becomes zero and the l.h.s. is zero only if the T, /i^ and don't depend on 
the (provided they don't depend on the momentum p^). The 1-st order space-time deriva- 
tives of the T, /ifc, in the l.h.s. should be cancelled by the first nonvanishing contribution 
in the r.h.s. This means that the should be proportional to the 1-st order space-time 
derivatives of the T, /i^, U^. The covariant time derivatives D can be expressed through 
the covariant spacial derivatives by means of approximate hydrodynamical equations, valid 
at the same order in the gradients expansion. Let's derive them. Integrating the (14.191) over 
the (^2n)^'p^ with the fk = fk^'' in the l.h.s. with inelastic collision integrals retained and using 
the (I4.24P and the (14. 3 p one would get (which can be justified with explicit form of inelastic 
collision integrals) 

9^Arf )^ = Dnk + rikV.U^ = h, (4.27) 

where Ik is the sum of inelastic collision integrals integrated over momentum. It is responsible 
for nonconservation of the total particle number of the k-th particle species and has property 
J2k1akh = 0. If (^^"^'[/fe] = then h = which results in conservation of the total particle 
numbers of each particle species. Multiplying the f l4.27p on the qak and summing over k one 
gets the continuity equations for the net charge flows: 

a^ATf/^ = Dn, + n.W^U^ = 0. (4.28) 
Then integrating the (I4.19p over the (2t)3pQ "^i^^ /fc — fk*^ l.h.s. one gets 

g^T^o)pu ^ Q^^^uPU"" - PAP") = 0. (4.29) 

There is zero in the r.h.s. even if inelastic collision integrals are retained because they respect 
energy conservation too. Note that the Boltzmann equations (I4.19P permit self-consistent 
consideration only of the ideal gas energy-momentum tensor and net charge flows. After 
convolution of the (I4.29P with the one gets the Euler's equation 

DU^ = -^—VP = —VP. (4.30) 

e + P hn ^ ^ 

After convolution of the (I4.29P with the Ui, one gets equation for the energy density 

De = -(e + P)V^JJ^ = -hnV^U^. (4.31) 

To proceed further one needs to expand the l.h.s. of the Boltzmann equations (I4.19P over 
the gradients of thermodjTiamic variables and the flow velocity. Let's choose /Xq and T as 
independent thermodynamic variables. Then for the Dfj!^^ one can write expansion 

Qf(o) ;^f(o) r)f(°) 

Writing the expansion for the Dua and the De one gets from the (I4.28P and the (14.310 : 

Dn^ = J2 ^Df^, + ^DT = -n^V.UP, (4.33) 
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De = ^DT + J2 = -hnV^U''. (4.34) 

a 

The solution to the system of equations (14.331) . (I4.34p can easily be found: 

DT = -RTV^U^", (4.35) 

D/x, = TY,Kb{RBh - X6)V^f/^ (4.36) 



where 



R ^ . (4.37) 



b 



and 

dua n ~ dua n de de 

^-^^a., ^-^5., ^--^M, dirr ^'-''^ 

Above it is assumed that the matrix Aab is not degenerate, which is the case usually. Other- 
wise uncertainties or singularities from -D/ia enter in the bulk viscosity. Using the ideal gas 
formulas (14. lip and (I4.12p one gets 



^ab = ''^qakqbkXk, Ea = '^^qakXkek, Ba = Eg — Ag, 



bfJ'b, 



k 



= ^Xk{'ihk + zl- fikh) = ^Xk{'ihk + zl) -^Egfig. (4.39) 

k k a 

From the (14.391) one can see that the matrix Ag^ has positive diagonal elements and in the case 
of one kind of charge it's always not degenerate. For the special case of vanishing chemical 
potentials, = 0, the quantities n^, Xg, Bg, Eg tend to zero because the contributions 
from particles and anti-particles cancel each other and chargeless particles don't contribute. 
Then from the (14.350 and the (14.360 one finds 

DT\,a=o = -^V^t/^ (4.40) 



Df^aU^=o = 0. (4.41) 

This means that for vanishing chemical potentials one can simply exclude them from the 
distribution functions (if one does not study diffusion and thermal conductivity). In systems 
with only elastic collisions each particle has its own charge so that one takes qgk = Sgk and 
gets 

Aki = SkiXk, Bk = Xk{ek — fik), Ek = CkXk, R = — 

a,b k k 



13 



Then the equation for the DT fl4.35p remains the same with a new R from the fl4.42p and 
the equations (I4.36P become: 



(4.43) 



Note that in systems with only elastic collisions the D^^ does not tend to zero for vanishing 
chemical potentials so that the /ifc could not be omitted in the distribution functions in this 
case. Because the heat conductivity and diffusion are not considered in this paper their 
nonequilibrium gradients are taken equal to zero, V uP = ^ uT = V^Ha = 0. Using the 
(05]) . K36\f and f OOj) the l.h.s. of the (KWf can be transformed as 

iPtU.D + py,)fi'^ = -T/f vr,^< + T/fg, V,[/^ (4.44) 

where 

a,b 



Qk = r^^-R]+Tk 



- ^4- (4.45) 



Using the (14.171) one can notice that the useful equality vr^vr^V^f/i, = Tr^Tr^V^f/;^ holds. In 
systems with only elastic collisions the Qk simplifies in agreement with [36|: 

= - 7) r', + rfc((7 - l)k - 7) - (4.46) 

where the assignments 7 from [36] is used. It can be expressed through the c„, defined in 
the (I4.42p . as 7 = ^ + 1. Introducing symmetric round brackets 

and assignments 



al ^ {Qk, inY), Ik = ((rfc)X<,^f^), < = (1, {ny)k, (4.48) 

and using explicit expressions of the from Appendix |A] one finds for the and the in 
systems with elastic and inelastic collisions: 

al = l + Y^ QakA-J {RBb - Xb) - [ik -fik)R, (4.49) 

a,b 

al = hk + Yl ekQakA-^iRBb - Xb) - {?>hk + zf + //^(l - hk))R. (4.50) 

a,b 

Then using the ( 14.49P and the ( ]4.50p one can show that 

'^qakXkal = 0, (4.51) 
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J2xkal = 0. (4.52) 

k 

o 

Because the gradients Vfj,U^ and V^U^ are independent the (I4.5ip and the fl4.52p are direct 
consequences of the local net charge fl4.28p and the energy-momentum fl4.29p conservations. 



Quantities (l,7r^7r^) and (Pfc,7r^7r^) vanish automatically because of the special tensorial 



structure of the vr^vr^."'^'^ 



The next step is to transform the r.h.s of the Boltzmann equations fl4.19p . After substi- 
tution of the (14.261) in the r.h.s. of the (I4.19p the collision integral becomes linear and one 

gets 

C^r[A]^-/rE^^'H' (4.53) 



where 



J^tM = 7^ / fS\^k + ^11 - A - ^'uWku (4.54) 



The unknown function ip^ is sought in the form 

1 



= (^-Afc(Pfc)V^t/^ + CkipkKir^k'^^U.j , (4.55) 

where (t{T) is some formal averaged cross section, used to come to dimensionless quantities. 



Then using the fl4.44p and the fl4.53p and the fact that the gradients V ^JJ^ and V^t/i, are 
independent the Boltzmann equations can be written as independent integral equations: 

Qu = Y.xiLti[Akl (4.56) 
I 

Y.xiLt[CuK<\. (4.57) 



I 

where the dimensionless collision integral is introduced 

Ltixk] = —^Xl\[Xk]. (4.58) 
niTa{T) 

In case of present inelastic processes the l.h.s of the f l4.56p is set by the source term f l4.45p 
and the r.h.s. contains linear inelastic collision integrals. After introduction of inelastic 
processes the source term in the f l4.56p becomes much larger as demonstrated in Sec. IVI A[ 
Using the equations f l4.36p and f l4.35p and the ideal gas formulas f l4.39p one can check that in 
the zero masses limit the Qk fl4.45p tend to zero and Dfia = that is the fia don't scale and 
the distribution functions become scale invariant. The source term of the shear viscosity in 
the f l4.57p doesn't depend on the presence of inelastic processes in the system and originates 
from the free propagation term ^^fe in the Boltzmann equation. 



o o 

Direct computation gives (1,^^)^ « {CiVUP + C2^''P)^''\p = 0, (p^tt^O^ oc [CiU^VUP + 
CiU^A'^P + dU" A^P)/^^"'„p = 0. 
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B. The transport coefficients and their properties 



After substitution of the /^^^ with the fl4.55p into (14.41) and comparison with the (I4.16P 
one finds the formula for the bulk viscosity 

1 T 

e = -3 ^ 5Z a;fc(A'^"7r^,fe7r,,fc, Ak)k, (4.59) 



k 



and for the shear viscosity 

lOa(r) 



I T sr^ ° o 

^ " -^—rnl^^k{T^kK^'^kT^kt,T^ku)k, (4.60) 



k 



where the relation A'^^^^A^A^ = 5 is used. 



In kinetics the conditions, that the nonequilibrium perturbations of distribution functions 
does not contribute to the net charge and the energy-momentum densities, are used as 
convenient choice and are called conditions of fit. They reproduce the Landau-Lifshitz 
condition (I4.15p . The conditions of fit for the net charge densities can be written as 



T.1-"'j0iflAU4'^'^>^ = a, (4,61) 

and for the energy-momentum density can be written as 

d^Pk 



o 

For the special tensorial functions CkTfkfJ^ in the (I4.55P they are satisfied automatically 
and for the scalar functions Ak they can be rewritten in the form (the 3-vector part of the 
fl4.62p is automatically satisfied) 

y^^qakXkiTk,Ak)k = 0, ^Xfc(r^, Afc)fc = 0. (4.63) 



For a single-component gas with only elastic processes the conditions (14.611) and fl4.62p 
exclude zero modes that is nonphysical solutions, proportional to the and a constant, for 
which elastic collision integral vanishes. In a single- component gas with inelastic collisions a 
constant is not a zero mode. For a multi-component gas these conditions would just modify 
the functional space on which solutions are sought. With help of these conditions of fit one 
can show explicitly essential positiveness of the ^. Namely, using the conditions of fit (I4.63p . 
the equation (I4.56P and the identity A'^'^vr^ fevr,^ ^ = ^1 ~ '^k ^^e bulk viscosity (I4.59P can be 
rewritten as 

where the square brackets are introduced for sets of equal lengths {F} = (Fi, F^, ...), 
{G} = (Gi, Gk, ■■■)■ 

{G}] ^ E 7^ / fffniFk + Fn - F\ - F'u)G,W,,. (4.65) 
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Using the time reversibility property of the Wm one can show that equahty 

{Fk + Fu - F'k - F'u)Gu = \{Fk + F^i - F\ - F\i){Gk + Gu - G\ - G\i), (4.66) 

holds under integration and summation in the (14. 65^ . Then one gets the direct consequence 

[{F}, {G}] = [{G}, {F}], [{F}, {F}] > 0. (4.67) 

This proves the essential positiveness of the ^. Similarly using the (I4.57P the shear viscosity 
can be rewritten in essentially positive form 

I T ^ f^— „ \ IT 



(^kK^CkT^kp^Tik^ ^ = \ y2^iLl\[Gk7rM^Cknk^,7iku 



10 (t(T) , 

^ ' k ^ ^ 1^ ^ ' k \ I / k 

= ^^[{^^}'{^^}]- (4.68) 

The considered variational method allows to find approximate solution of the integral 
equations (14.56^ and (I4.57P in the form of linear combination of test-functions with coeffi- 
cients, found from condition to deliver extremum to some functional, which is equivalent to 
solving the integral equations. One could take this functional in the form of some special 
norm as in 361. Or one can take somewhat different functional like in 41 , which is more 



convenient, and get the same result. This generalized functional can be written in the form 
i^M = E^'=(^^'^^'^^-)^-^[{^^''}'i^'^M...4], (4.69) 



where S'^'"'^ = Qk and XkiJ....v = Ak for the bulk viscosity and S^'"" = vr^vr^, Xk'"" = ^k^^k'^k 
for the shear viscosity. Equating to zero the first variance of the (I4.69P over the Xkii...u one 
gets 

Yl ^^(Sr, Sxk,...u)k - [{xr}, {5Xk,...u}] = 0. (4.70) 

k 

Because variations 6xkfi...u are arbitrary and independent the generalized integral equation 
follows then: 

Sr = T.''iLl\[xk,....]. (4.71) 

I 

The second variation of the f l4.69p is 

S'F[x] = -[{Sxr}ASXk,....}]<0, (4.72) 

which means that solution of the integral equations fl4.56p and f l4.57p is reduced to variational 
problem of finding the maximum of the functional ( I4.69p . Using the (I4.7ip the maximal value 
of the (14.690 can be written as 

Fmaxlx] = ^[{Xk"''}, {XfcM...J]lx=Xmax- (4-73) 
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Then using the (14.641) and the fl4.68p one can write the bulk and the shear viscosities through 
the maximal value of the F[x\ 

^ ~ , (4.74) 

--Ak 



-F 

' /m\ max 

oiT) 



—QkiXkfi...iy—^k 



1 T 



V 



5a(T)- 



(4.75) 



This means that the precise solution of the fl4.7ip delivers the maximal values for the trans- 
port coefficients. 

The approximate solution of the system of the integral equations fl4.56p and f l4.57p are 
sought in the form 

A, = Y,AlTi:, (4.76) 

r=0 



n2 



r r 



(4.77) 



where ni and n2 set the number of used test-functions. Test-functions used in [4l| would 
cause less significant digit cancellation in numerical calculations but there is a need to reduce 
the dimension of the 12-dimensional integrals from these test-functions as more as possible 
to perform calculations in reasonable time. The test-functions in the form of just powers of 
momentums seem to be the most convenient for this purpose. Questions concerning unique- 
ness and existence of solution and convergence of the approximate solution to the precise one 



are covered in 36 



As long as particles of the same particle species but with different spin 
states are undistinguishable their functions fl4.55p are equal and the variational problem 
is reduced to variation of coefficients Al, and CI, and the bulk fl4.64p and the shear fl4.68p 



viscosities can be rewritten as 



T 



N' 



"1 



a(T) ^ ^ 

^ ' k'=l r=0 



(4.78) 



1 T 
Waif) 



N' 712 



(4.79) 



=1 r=0 



After substitution of the approximate functions A^, f l4.76p and C^/ fl4.77p into the fl4.69p and 
equating the first variation of the functional to zero one gets the following matrix equations 
(with multi-indexes (/', s) and {k', r)) for the bulk and the shear viscosities correspondingly^'^ 



N' m 



(4.80) 



V=l s=0 



One can first derive the same equations for the Ak and Cfe, treating them as different functions for all fc, 
with the coefficients and CJf having the same form as the A^^i, and C|!f;,. Then after summation of 
equations over spin states of identical particles and taking — Ak' , Ck = Ck' one reproduces the system 
of equations for the Ak' and Ck' ■ 
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N' 712 

^i^'ll' = Y.Y.^i'^'^v, (4.81) 

Z'=l s=0 

where introduced coefficients A]^i, and C^f;/ are 

Al%, = Xk'Xv[T\ Tl]kn' + ^k'l'Xk' ^ Xm'ir'', T%irn', (4.82) 

m' 

° o _ — ^ ° o 

m' 

They are expressed through the colhsion brackets 

^i]'^' = T6a ^2 2 2 ^ / e--''--" (F, - F',)Gi,iyH. (4.84) 

T\AnYzizfK2[zk)K2{zi)a[T) 

The colhsion brackets [F, G]fci are obtained from the last formula by replacement of the Gu 
on the Gk- Due to time reversibility property of the Wki one can replace the Gu on the 
1(^1^ — G'li) in the (14.841) . Then one can see that 

[t\t%i>Q. (4.85) 

Also it's easy to notice the following symmetries 

[F, Gi]ki = [G, F^]ik, [F, G]ki = [G, F]^. (4.86) 

They result in the following symmetric properties A^f^, = Af/^,, G^f^, = Gf,^/. Also the micro- 
scopical particle number and energy conservation laws imply for the ^4^4'- 

Al% = 0, (4.87) 
J2a11, = 0. (4.88) 



The (14.871) together with the = (1A13P means that the equations with r = in the (I4.80p 



are excluded. From the (14.881) it follows that each one equation with r = 1 in the (I4.80p 
can be expressed through the sum of the other ones, reducing the rank of the matrix on 1. 
To solve the matrix equation (14.801) one eliminates one equation, for example with k' = 1, 
r = 1. One of coefficients of Aj, is independent, for example, let it be A\,. Using the (I4.88p 
the matrix equation (I4.80p can be rewritten as 

TV' N' m 

l'=2 l' = l s=2 

Then using the (lAlSp and the (14.521) the bulk viscosity (I4.78P becomes 



J. N' J, N' rii 



^ ' k'=2 ^ ' k'=l r=2 
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Then the coefficient A\, can be ehminated by shift of other Aj, and be imphcitly used to 
satisfy one energy conservation condition of fit. The particle number conservation conditions 
of fit are imphcitly satisfied by means of the coefficients A^,. The ffist term in the fl4.90p is 
present only in mixtures. That's why it is small in gases with close to each other masses of 
particles of different species (like in the pion gas). In gases with very different masses (like 
in the hadron gas) contribution of the ffist term in the (I4.90p can become dominant. 

Analytical expressions for some lowest orders collision brackets which enter in the matrix 
equations f l4.8ip and fl4.89p can be found in Appendix Higher orders are not presented 
because of theirs bulky form. 



V. THE NUMERICAL CALCULATIONS 

The numerical calculations for the hadron gas involve roughly 2 ^'^^^ 12-dimensional 
integrals, where N' is the number of particle species and n is the number of used test- 
functions (called the order of calculations). The 12-dimensional integrals being the collision 
brackets [F, Gi]ki and [F, G]ki can be reduced to 1-dimensional integrals, expressible through 
special functions. To compute these special functions precisely and fast they were expanded 
into series at several points. This allows to perform calculations at the 3-rd and the 6-th 
orders for the shear and the bulk viscosities correspondingly. Because analytical expressions 
for the collision brackets are bulky the Mathematica j42| was used for symbolical and some 
numerical manipulations. The numerical calculations are done also for temperatures above 
T = 140 MeV, where inelastic collisions are expected to play important role, for future 
comparisons and for the bulk viscosity to show the position of its maximum in the hadron 
gas. 



The hadron-resonance particle list from j43| is used in calculations at zero chemical po 



tentials. The results for the shear and the bulk viscosities are shown in fig. |2] and fig. [3] 
correspondingly. The maximal relative errors are less than 0.172% for the bulk viscosity and 
less than 0.444% for the shear viscosity. As one can see in the fig. [T]the bulk viscosity of 
the pion gas has maximum approximately at the temperature equal to the half of the pion's 
mass. In the hadron gas this maximum shifts towards the value T ^ 200 MeV. 

The hadron gas mass spectrum dependence of transport coefficients is investigated. Par- 
ticles in the list were sorted over their masses and the list was cut on the 3 GeV, 2 GeV, 
1 GeV and on the pion's mass. It was found that mass cut can be done on the 3 GeV 
which results in negligible errors. The particle list, cut on the 3 GeV, contains 462 particle 
species with anti-particles included. The comparison of the shear viscosities calculated with 
different mass cuts is depicted in fig. |H As one can see the shear viscosity changes no 
more than in 1.5 times. And the bulk viscosity mass spectrum dependence is very strong as 
depicted in fig. O The ratio of the bulk viscosity of the hadron gas to the bulk viscosity of 
the pion gas reaches 8 at T = 120 MeV and 13.7 at T = 140 MeV. 

The ratio of the shear viscosity to the entropy density rj/s and the ratio of the bulk 
viscosity to the entropy density ^/s in the hadron gas is shown in fig. [61 The ratio ^/s 
doesn't have a maximum and is descending function of the temperature. The entropy 
density is calculated by the formula ( \B7\i using the ideal gas formulas in the (14. lip and the 
f l4.14p . The ratio of the bulk viscosity to the shear viscosity is shown in fig. [71 

The dependence of the rj/s and the ^/s from the temperature, calculated along the freeze- 
out line, is found too and is depicted in fig. [SI As was discussed in Sec. [TTH calculations 
with large chemical potential may contain large deviations especially for the bulk viscosity 
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FIG. 2: The shear viscosity as function of the temperature in the hadron gas calculated up to the 
3-rd order. Chemical potentials are equal to zero. 
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FIG. 3: The bulk viscosity as function of the temperature in the hadron gas calculated up to the 
6-th order. Chemical potentials are equal to zero. 

and are rather estimating. At considered collision energies strange particle numbers are not 
described well with the statistical model. It's expected that this is because they doesn't reach 
chemical equilibrium before the chemical freeze-out takes place. After introduction of strange 
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FIG. 4: The ratio of the shear viscosity Tj^ceV with the mass spectrum, cut on the 3 GeV, to the 
shear viscosities with the mass spectrum, cut on the 2 GeV, 1 GeV and on the pion's mass, as 
function of the temperature. Calculations are done for zero chemical potentials. 
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FIG. 5: The ratio of the bulk viscosity S,3GeV with the mass spectrum, cut on the 3 GeV, to the 
bulk viscosities with the mass spectrum, cut on the 2 GeV, 1 GeV and on the pion's mass, as 
function of the temperature. Calculations are done for zero chemical potentials. 



22 




0.08 



0.12 0.16 

T (GeV) 



0.20 



0.24 



FIG. 6: The ratio of the shear viscosity to the entropy density and the ratio of the bulk viscosity 
times 100 to the entropy density as functions of the temperature in the hadron gas. Calculations 
are done at zero chemical potentials. 
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FIG. 7: The ratio of the bulk viscosity to the shear viscosity times 100 as function of the temper- 
ature in the hadron gas. Calculations are done at zero chemical potentials. 



saturation factors 7^ experimental data gets described well [37[, [38[. These calculations 
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were done using the old particle list from the THERMITS package |4J] (without charmed 
and bottomed particles, which comprise 358 particle species including anti-particles). All 
variables' values of freeze-out line including the strangeness saturation factor 7^ see in (23| . 
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FIG. 8: The shear viscosity and the bulk viscosity times 100 as functions of the freeze-out values 
of the temperature T. All other freeze-out line variables' values can be found in 23l |. 



VI. ANALYTICAL RESULTS 
A. The single-component gas 

In the single- component gas with one test-function the matrix equations can be easily 
solved and the shear (14.791) and the bulk (14.901) viscosities become (indexes " 1" of the particle 
species are omitted) 

' (6.1) 



10 (t(T) coo 



In this approximation the explicit closed-form (expressed through special and elementary 
functions) relativistic formulas for the bulk and the shear viscosities were obtained in 26 



for the model with constant cross section with the ideal gas equation of state. There the 
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parameter a = 2r. In [36[ they are written through the parameter a = 2r . The results 



are: 



^ 647r r2 {5z^ + 2)K2{2z) + {3z^ + 49z)K3{2z) ' ^ ' ' 



1 T z'Kl{z)[{5-3^)h-3^]' 
^ QAnr^ 2K2{2z) + zKs{2z) ' ^ ' ' 

where 7 = — + 1 = f'^^^~!!'^ . Though the correct resuh for the shear viscosity is 
^ 647rr'^{15z^ + 2)K2{2z) + {3z^ + A9z)K3{2z)' 



This result is in agreement with the result in [15|, |45l] . To get the fl6.5p and f l6.4p the collision 
brackets in the fl4.83p and A^"^ fl4.82p can be taken from Appendix O with Zk = zi = z 
and the 7° and can be taken from Appendix \M In the nonrelativistic limit z ^ 1 one 

gets 

5 T ^/ 25 1 \ , , 

25 T / 183 1 \ ,^ 
In the ultrarelativistic limit 2; ^ 1 one gets 



^ ^ 2^?^^' + (i " """^ + ''^'') + + . (6.9) 

where 'Je is the Euler's constant, •je ~ 0.577. 

The perturbation of the distribution function cp f l4.55p can be found too: 

<^ = -^p^ ( + AW + AW^)V^W' + C°^F^V74^ , (6.10) 



na(r) 

where the C° is equal to 



,^l^aiT)_ z'Kl{z)h 

647r r2 (15^2 + 2)/s:2(2-2) + (3^3 + 49^)7^3(2^) ' ^ ' ' 



It is the differential cross section for identical particles. The total cross section is / ^2r^ = Anr^ 
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and the is equal to 



1 cjiT)z^Kl{z)[{5-3^)h-3j] 



647r r2 2^2(2^) + ^^^3(2^) ' ^ ' ^ 



The and the are used to satisfy the conditions of fit fl4.63p and are equal to 



where the a* can be found in Appendix |Al In the nonrelativistic limit one has 

= ^rr.. o o -(^' + - ^')V^f/^ + 27r/^7r-V,.f/. . (6.14) 



32y2T' 



In the ultrarelativistic limit z <^ 1 one has 



^ = 4^^^ (^-5^2(r2 + 8r - 12)V^f/^ + 36n^n^V j . (6.15) 

Note that although the shear viscosity diverges for T — )■ 00 the perturbative expansion over 
the gradients does not break down because the (p does not diverge and tends to zero con- 
versely. 

The phenomenological formula, coming from the momentum transfer considerations in 
the kinetic-molecular theory, for the shear viscosity is rjph oc ln{\p\) (with the coefficient of 
proportionality of order 1), where is the average relativistic momentum (13.31) . / is the 
mean free path. It gives correct leading m and T parameter dependence of the (16. 5p with 
quite precise coefficient. The mean free path can be estimated as / ~ l/{(rtotn) (see Appendix 
[pi). Choosing the coefficient of proportionality to match the nonrelativistic limit one gets 



5 V^K,/2{m/T) 
''^''"647^ r2 K2im/T) ' ^ ' 

If the bulk viscosity is expressed as oc the coefficient of proportionality is not of 

order 1. In the nonrelativistic limit it is 25/(512a/22;2) and in the ultrarelativistic limit it 
is ^''/(8647r). To reproduce these asymptotical dependencies the bulk viscosity should be 
proportional to the second power of the averaged product of the source term Q and the r 
that is to the (a^)^. 

If the system has no charges then terms proportional to the Tk in the (I4.45P are absent. 
This results in very different values of a^. In particular for the single-component gas in the 
case z ^ 1 one gets 



Q!2|gil=l 15 



+ .... (6.17) 



and in the case z ^ 1 one gets 



gii=o 



13 + .... (6.18) 

l<7ii=l 

In the case z ^ 1 this means for the bulk viscosity that it becomes divergent over the 
temperature at least as 0{Tz^^^'^) oc 0{T~^^^'^). However this dependence holds if number- 
changing processes are not negligible. 
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B. The binary mixture 

The mixture of two species with masses mi, 1712 with different classical elastic differential 
constant cross sections af^, crfa = c^i, C22 is considered in this section. Using the fl4.79p with 
n2 = and solving the matrix equation fl4.8ip one has for the shear viscosity: 

V = T5^i;[(^i'^i°)'^°'°' - 2xvX2.7?72°Ci% + (^2'72°)'C',%], (6.19) 

where = C^?yC^?2' - {Cv2'f- The collision brackets for the C^i], (^M) can be found in 
Appendix O and the 7° can be found in Appendix |A1 

In important limiting case when one mass is large ^2 3> 1 {g2 and jl2 are finite so that 
X2' <^ 1) and another mass is finite one can perform asymptotic expansion of special func- 
tions. Then one has xi' oc 0(1), 7? oc 0(1), X2' oc 0{e~^^Z2^'^), 7° oc 0(^2)- Collisions of light 
and heavy particles dominate over collisions of heavy and heavy particles in the O2/2' and one 

has [vr/^vr^, 7f^^]2i oc 0(2:2), O^/^ / oc 0{e~^^Z2 )■ In the Oj'^, collisions of light and light par- 

° o 3 ,2 

tides dominate and one gets oc 0(1). And [^F^, vfl^Ivfl^] 12 oc 0(1), 0°,°' oc 0{e-'^Z2'"). 
In the shear viscosity the first nonvanishing contribution is the single-component shear vis- 
cosity (16. 5p . where one should take = al[ and z = zi. The next correction is 

Z1-A1+A2 



64:V2n{3 + 3zi + zl)gia 



The approximate formula 23l 



V = y^^VkXk, (6.21) 

k 

where rj^ is given by the (16. 5p or the (I6.16P with mass and cross section al\, would give 
somewhat different heavy mass power dependence 0(e~^2 2;|). 

Using the ( I4.90p with rii = 1 and solving the matrix equation ( I4.89P one has for the bulk 
viscosity: 

T (a;2'«2)^ T xiiX2ia\ci\ 

Using definition of the Al}2f (14.821) and the fact [T,Ti]ki + [T,T]ki = (lC32p one gets 
Al}2, = xi'X2'[t,t]i2. Using the (14.851) one gets [t, r]i2 > 0. Then using xi'a\ + X2ia\ = 0, 
coming from the (14.521) . the bulk viscosity can be rewritten as 

^ = ^7^ — \ ^ = — \ ^ > 0- (6.23) 

a{T) xr[T,T\i2 a{T) X2'[t,t\i2 

The collision bracket [T,r] 12 can be found in Appendix O and the al can be found in 
Appendix [XI 

In the limiting case ^2 3> 1 one has Xi' oc 0(1), X2' oc 0(e~^^Z2^^)) '^1 0^ 0(e"^^2;2^^), 
al oc 0(1), y422 oc All 0{e~^^ zl^'^), [r, r]i2 oc 0{z2^). Then for the bulk viscosity one gets 

. 5/2 g2Te'^'^+^'^+'^2z'^ -5-2h'^, +10hif , 

P = e-'^zl' — ^ . ^ . — + .... 6.24 

Vm^gxaf2{zl^?>zi^'i)\z\-\-h\^^hif 
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VII. CONCLUDING REMARKS 



The shear and the bulk viscosities of the hadron gas and the pion gas were calculated at 
low temperatures in the model with constant cross sections. The physics of the bulk viscosity 
is very interesting. In particular it was found that in mixtures with only elastic collisions it 
can strongly depend on the mass spectrum. For instance, at temperature T = 120 MeV the 
bulk viscosity of the hadron gas is larger in 8 times than the bulk viscosity of the pion gas. 
Also the bulk viscosity can strongly depend on quantum statistics corrections, equation of 
state and inelastic processes which can be explained by nontrivial form of its source term. 
Inclusion of inelastic processes in pion gas at T = 140 MeV results in increase of the bulk 



viscosity roughly in 44 times according to comparison with results of the paper [20|. It's 
a future task to switch off carefully inelastic processes where they can be considered as 
negligible ones to perform calculation of the bulk viscosity in the pion gas and the hadron 
gas. The shear viscosity is less dependent on the mass spectrum and on quantum statistics 
corrections and its source term is some trivial function which doesn't depend on inelastic 
processes. 
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Appendix A: Values of a^, jj, and 
Their definitions are 



al = {Qk,{rkY), Yk = {{rkyKK,rf^), 4 = {1, {rky)k, (Al) 
where the round brackets are 

(F,G)fc^ . oj. / F{pk)G{p,)e-^-. (A2) 



Pk 



Then one can rewrite the a|. as 



There is recurrence relation for the a|: 



rdr{r'-zlf''r^e-\ (A3) 



al = {s + l)al-' + zlat' - {s - 2)zlal-\ (A4) 
It can be derived from the flA3D written in the form 



zlK^izk) A, ' ' ' 
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Then after integration by parts the recurrence relation follow. Some values of the are 



2^ 


(A6) 


4 = 1, 


(A7) 


(^k = hk — 1, 


(A8) 


al = 3hk + zl, 


(A9) 


4 = (15 + zDh + 2zl 


(AlO) 


al = 6{15 + zl)h + zl{15 + zl), 


(All) 


= (630 + A5zl + zl)hk + 5^^(21 + zl). 


(A12) 



The can be expressed through the a'^ after integration of the (14.451) (or the fl4.46p if only 
elastic collisions are considered) over momentum and using the definition flAip . For systems 
with only elastic collisions some values of the are written below, in agreement with [36i] : 

al = 0, (A13) 

1 2(c.-9)/ifc + 3/ig-32:g 7fc - 7 

«fc = = 7, (A14) 

Cv 7fe - 1 

al = 2h - 3^ - 3^^ = (5 - 37)/.. - 37.^^, (A15) 

where the assignments 7 and 7. from 36|] are used. They can be expressed through the c„ 
and the c„,fc, defined in the (14. 42 p . as 

7 = - + 1, lk = — + l- (A16) 



The 7^ can be rewritten as 



SzlK^izk] 
Then it can be rewritten through the a 



2 1 f°° 

= oT^^TTTT / dr{r' - ziy/'r^e-^ (A17) 



Some values of the 7^ are 



Yk = l{a^^'-2zla^+' + zlan. (A18) 



7° = lOh, (A19) 
yl = 10{6hk + zl), (A20) 
yl = 10{7zl + h{42 + zl)). (A21) 
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Appendix B: The entropy density formula 

The Gibbs's potential is defined as 

^{P,T) = E{S,V) - ST + PV. (Bl) 
The differential of energy is defined as 

dE = TdS - PdV + l^kdNk = TdS - PdV + f^adNa, (B2) 

k a 

where it is rewritten through the independent chemical potentials and the particle net 
charges Na. The differential of the $ then reads: 

d$ = -SdT + VdP + J2 f^adNa. (B3) 

a 

Because the $ is function of intrinsic variables P, T and extrinsic Na the only possible form 
of it in the thermodynamic limit is 

$ = 5^iV,0,(P,T), (B4) 

a 

where 0a are unknown functions. Then from the ( ]B3I) one gets = fj^a, which means that 
4'a = l^a- Then substituting the (lB4p into the (IBip one gets the relation 

NafiaiP, T) = E{S, V)-ST + PV. (B5) 

a 

Being written for local infinitesimal volume it transforms into the expression 

J^n./i, = e-sT + P, (B6) 



from where the entropy density s can be found: 

e + P 



T 

a 



ngftg. (B7) 



Appendix C: Calculation of the collision brackets 

The momentum parametrization and the most transformations of the 12-dimensional inte- 
grals used below are taken from 36[ (chap. XI and XIII). Let's start from some assignments. 



The full momentum is 

p^^p',+p';,=p",+ru- (ci) 

The "relative" momentums before collision Q^^ and after collision Q'^ are defined as 

= A^'^(p,, - pu.), Q'^ = A'i^'ip',, - p\i,), (C2) 



30 



with the assignment 



pupu 
p2 ' 



(C3) 



where = P^P/j,. The covariant cosine of the scattering angle can be expressed through 
the Qf" and the Q"" as 



cos© = 



Q-Q' 



-QW-Q" 



(C4) 



where • denotes convolution of 4- vectors. One also has — Q' and 



= 4ml - (1 + c^kirP' = -{P'- Ml,) 

where 

Mki = mk + mi, Hki 



Mil 



kl 



Ml 



kl 



(C5) 



mkUii 



ml — m 



mk + mi' 



Oikl 



p2 



- = sign(mfe - mi)Jl 



Mm P' 



(C6) 



The function sign(a;) is equal to 1, if x > and equal to —1, if a; < 0. Note that not all 
and are independent: 



p^g^ = o, P''Q'^ = o. 



(C7) 



To come from the variables (p^jPii) to the variables {P^,Q^) in the measure of integration 
first one has to come from the (p^jPii) to the (p^ + Pii.p^ — Pii) (the determinant is equal to 
16) and then shift the relative momentum p^—p^^ on the akiP^^. Analogically for the {p'kiP'ii) 
and the (P'*, Q'^). The inverse relations for the Pk-,Pii-,p'k->P'ii through the P'^, Q^, Q'^ are 



Pk 



l{l + aki)P^ + lQ^, 



(C8) 



(C9) 



Pk = l{^ + C^kl)P'+lQ"': 



(CIO) 



There is a need to calculate the following integrals 

r{a,b,d,e,f\g,r) ^kl 



(Cll) 



J, 



kl 



6-^-^/^(1 + akiy 



2^2 



Wki. 



(C12) 
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After nontrivial transformations described in more details in 36| one arrives at (the constant 
cross section approximation is used) 

j{a,b,d,ej\q,r) _ TT[a + 6 + ijUCTikl ST \^ \^ \^ \^ f ^ M ^ ^2('^l+r■l+fc2+fc3) 



\Zk + ZiJ \2hJ \qij \rij \ k2 / 

/d±e 

xM ^^^^ J/f2(a + /-gi-ri-A;2-A;3) + 3,6 + ^-/i + l,Zfe + zz 
where 

min((i,e) 

^2"'^ = ^ E ^^''^^^(^^e,,), (C14) 

^ ' 9=0 

where al\ = ^kio'ki is the classical elastic differential constant cross section. The 
equal to nonzero value 

if the difference f — g is even and g < f. The e, (7) is equal to nonzero value 

= (rf-9)!!(<i + 9+l)!!(e-9)!.(e + s + l)r f^^'" 

if < min((i, e) and both c? — (7 and e — (yf are even (which also implies that c? + e is even). 
The [...] denotes the integer part. The integral / is 

/oo 
duu'^''Kn{xu). (C17) 

Also there is the following frequently used combination of the J integrals: 

j,{a,b,d,e,f\q,r) _ -^^u f f\ j^2^^^2(/-'Lt) j(a+fc,6,d+e,0,0|i3+2?i,r) _ j{a,b,d,ej\q,r) (C18) 

The first term in the difference is obtained by the replacement of the Q' on the Q. Using 
this fact the J' can be rewritten in the form 

f, , nil W^./) 5 ' [f>/2] 

k I z\ tcj z\ I) gj=Ori=0 fc2=0 A:3=0 h=0 

^ (i^f )"^"^^'' - Gi) (:) (;) t: ^ o <-^) 

xf 2 jj/2(a + /-gi-ri-fc2-fc3) + 3,6 + ^-/i + l,Zfc + zJ , 
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where 



id,ej) 
kl 



a 



a{T) 



mm((i,e) 

{d + e + 1) ^ K{d + e,0,0) - ^ K{d,e,g)a 

9=0 



U,9) 



(C20) 



There is recurrence relation for the integral / (1C17P |36l . |46 



J(r,n,x) = (r - l)(r + 2n - l)/(r - 2, n, x) + (r - + (C21) 

For calculations one needs only integrals /(r, n, x) with positive values of the n and odd 
values of the r. If r > — 2n + 1 the / integrals can be expressed through the Bessel functions 
Kn{x) using the ( ]C21|) , when r = 1 or r = —2n + 1. Then using the following recurrence 
relation for the Kn{x) j46j 



2n 

X 



(C22) 



the final result can be expressed through a couple of Bessel functions. If r < — 2n — 1 
then the recurrence relation (1C21I) blows if one tries to express the /(r, n, x) through the 
/(— 2n + l, n, x). Using the (10211) the / integrals with r < —2n — 1 can be expressed through 
the integrals G{n, x) 



G{n, x) = I{—2n — 1, n, x) 
There is recurrence relation for the G{n,x): 



X 



duu " Kn{xu). 



(C23) 



Gin.x) = iG{n-l,x)-x "^KJx)). 

2n 



(C24) 



It can be easily proved by integration by parts of the flC23P and using the following relation 
for the K„(x) m 



— Kn{x) = Kn{x) - Kn-l{x). 

OX X 



(C25) 



It is found that collision integrals have the simplest form if they are expressed through 
G{n, x) with n = 3 or n = 2 and the Bessel functions K^lx) and K2{x) or K2{x) and Ki{x). 
It was chosen to take G{x) = G{3, x) and K^i^x), K2{x). The G{x) can be expressed through 
the Meijer function [47] 



G{x) 
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1 

-3,0,0 



The needed scalar collision brackets can be expressed through the J' as 



/(0,r+s— M— ?j,M,D,0|r— it.s— sj) 
kl ' ■> 



(C26) 



(C27) 



1 

2r'+s / ^ / ^ 



r 



\J' 



/(O/r+s— M— f ,M,ti,0|r+s— D,0) 
kl ■> 



(C28) 
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and the needed tensorial collision brackets can be expressed as 



r^r ^u" t/ n-s ^ °„ — 1 i ^\ i ^ \ ( ^\s-n1( T/(2,ni+n2,r-m,s-n2,0|2+ni,2+n2) 

ni=0ra2=0 \ / \ / 

I 9 THl,ni+n2,r-ni,s-n2,l|l+"i,l+"2) , jl{0,ni+n2,r-n\,s-n2,2\n\,n2)\ 
^ -^-^ kl kl ) 

1 /^^ ^ "'^ V ^\s+^.-n2( T/(l,ni+n2,r+l-ni,s+l-n2,0|l+m,l+n2) 

ni=0n2=0 ^ ' ^ ' 

^ _^.(0,n,+„2,.+l-n„.+l-n2,l|n„n2)) ^ 2 ^,^2]^^ ^ 1^2[^.+2^ ^.]^^ 

+ ^^^[T^r^^]./-^4^^^^^^]fc^ (C29) 

[rW,r^^], = ^ E E (; ) (; ) (J'g— -'01^— ^-0) 

ni=0n2=0 V 2/ V 1/ 
— 2 J''''''''*^"^"^'''"'"^'*^"^'^'^"''"^"'""^''^^ + j/(0,ni+n2,r-ni,s-n2,2|ni+n2,0)\ 
r+1 s+1 / , t\ / , t\ 

/(l,ni+n2,r+l-ni,s+l-r!,2,0|2+ni+r!.2,0) 



2r+s+3 

ni=0 n2=0 

T/(0,ni+n2,r+l-m,s+l-n2,l|ni+n2,0)N , ^ r r+2 ^5+2] , ^ ^2r r+2 
kl ' '^^ ikl^'^^kV 1^ \kl 

+ \4[r'.r^''\i-\4[r',r%i. (C30) 

Below some lowest orders collision brackets are presented with the following notations (one 
constant cross section approximation is used below and a'j^i are taken equal to a{T)): 

~ _ K^^Zk + zi) ~ _ K2{zk + zi) ~ _ G{zk'^zi) 
-f^i = '77~i — \~rr~i — -f^2 = -77-, — 777-7 — -f^3 = 



i^2(^fc)i^2(^/)' K2{Zk)K2{Zi)' K2{Zk)K2{Zi)' 

Zki = Zk + zi, Zki = Zk — zi- (C31) 



For the scalar collision brackets one has: 

'-k'^l ^kl 

where 

j(i>i) _ 



-[r, n],, = [r, r],, = ,7;2;2^(PS''^i^i + PS''^i^2 + Pis^'^^i^s), (C32) 



and 



P,V^'^^ = -2Zu{zti + 44Z,^, - 2Z,^J, (C33) 
n?''^ = 4(3^fe', + 8) + ?>2zl,Zl + 8Z,^,, (C34) 
Pis''^ = -34^.'^ (C35) 



[r, r,2],, = [r^ n],, = T;^(Pi;f ^^^i + P^l?K, + P^^^i^^s), (C36) 

^ki 
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where 



and 



where 



and 



where 



and 



where 



Piif ) = 2Z,i{zliZki + 84 + lQzl,Zki + ^2zlZl, + IQzuZl - AOZ^i), (C37) 

P^li^ = -4iZki{Zl + 8) - Szl,{Zl + 8) - 164Z,,(^,', + 8) 

+ 164^2^(^2^ - 16) + ^z,iZl,{Zl - 16) - %Zti{Zl + 8), (C38) 

P'^.? = 4iZli, (C39) 



[r, r^]., = [r^ r],, = T;2j^(i^lf ^^^i + P^'^^K, + P^H^ K,l (C40) 
"^^fe^/ ^ki 



P^l? = 2Z,,(4;^fc; - 84 + 164;^;^/ - 324;^.'; + ^^ZkiZli + 40Z,^,), (C41) 

^^22'^ = -zlMZl + S) + ^zli{Zl + ^)-lQzlMZl + S) 

- iQ4iZUZli - 16) + Sz,iZli{Zli - 16) + 8Z',i{Zli + 8), (C42) 

P^l? = (C43) 



[r', r',U = ^^^^iPsi?K, + P^l^^K, + Plli'K.l (C44) 



P^l^^ = -2Z,,[zUZl, + 2) + 6zUllZli-32)- 72zlzUzl, + 8) 

+ 24ZUZI + 96)], (C45) 

pS? = 4iiZti + lOZl^ + 16) - 64(Z,^, - 56^2, + 256) 

+ lUzl^ZUbZl, - 32) - 48ZUl^Zl, + 32), (C46) 

PSi^ = -4iZli[zUzl - 6) - 6Zl,], (C47) 



Pi^Y^ = -2Zki[zUzli + 2) - 364Zfe, + 184(^2^ + 16) + 96zl,Zki{Zli - 10) 

+ 2izlZl{Zl + 56) - ^Sz,iZUZli + 20) - 24^,^^,', + 100)], (C49) 
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pS? = 4i{Zti + iOZl + 16) + 12z'MZl - 24) - Gz^Zti - 72Zl - 384) 
- imzlMZl + 40) - 484^2 (13^2 - 224) + QQzkiZl{lZl - 80) 
+ 48Z,Vl3Z,^, + 48), (C50) 

P'&? = -4iZli[zUZl - 6) + l2z,iZ,i - 6Zl]. (C51) 
And for the tensor colhsion brackets one has: 

\:^,W^U = (PtuK, + P^TK2 + P^K,), (C52) 

'"^h^i ^ki 



where 



and 



Pi°f = -2Zm[zU5ZI - 8) + 244(^2^ - 16) - lUz^ZliiZl^ + 8) 

+ ^SZti{Zli + 72)], (C53) 

pS = 4i{^Zt, - AOZl, - 64) - 2AzU5Zti + 8Zl, + 128) 

+ ^^^4iZli{Zli - 16) - lQ2Zl,{hZl + 16), (C54) 

= -^4iZli[zUZl - 24) - 2AZU, (C55) 



W^. = ^irA^l^Sf ^1 + P^2K2 + pSKs): (C56) 

'^^fc^; ^ki 



where 



Pg? = 2Z,,[4(8 - 5Zli) + 72zU^Zl - 8) - 4804Z,K^,', - 4) 

- SS6zl,ZUzli + 8) + 2A0zmZUzI + 8) + 192Z,^,(Z,2, + 67)], (C57) 

P^T = 4(5^fc^ - ^OZl^ - 64) + 2404zf, - 244(5Z,^, + 48^^^ - 192) 

+ 19204Z,K^fe'; - 8) - 1924^2 (17^2^ - 112) + 1920zkiZli{Zi, - 8) 

+ 768ZU5Zl, + 6), (C58) 

423^ = -HiZliizUZl - 24) + A8z,iZ,, - 2AZl]. (C59) 

If Zk = zi then the G{x) function is ehminated everywhere and colhsion brackets simplify 
considerably. 
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Appendix D: Collision rates and mean free paths 



The quantity o o^'°/o ,o d'^p' yd^p'uh which enters in the elastic colhsion integral fl4.22p . 

Pfc'^li'P k'P w 

represents the probability of scattering per unit time times unit volume for two par- 
ticles, which had momentums |4' and pu' before scattering and momentums in ranges 
{p' k'lP' k' + dp' and (p' n' , p' n' + dp\i,) after scattering. The quantity j^^fk' represents 
the number of particles per unit volume with momentums in the range {pk',Pk' + dpk'). The 
number of collisions of particles of the k'-th species with particles of the Z'-th species per 
unit time per unit volume is then^^ 

Dei _^ ^ ll'l' [ d^Pk' d^pw d^p'^, d^p\.(^) 
l-^TTJ J p^, p-^i, p j^, p 

To get the corresponding number of collisions of particles of the k'-th species with particles 
of the I'-th species per unit time per particle of the k'-th species, Rf^^, one has to divide 
the fIDll) on the 'yk'i'fik' (recall that n^/ oc gk' by definition), which is the number of particles 
of the k'-th species per unit volume divided on the number of particles of the fc'-th species 
taking part in the given type of reaction (2 for binary elastic collisions, if particles are 
identical, and 1 otherwise). This rate can be directly obtained averaging the collision rate 
with fixed momentum pk of the k-th particle species 

Tyel — f d^Pll' ,3 , ,3 , „(0) Wkl' , . 

l^^J PkPll'P kP 11' 



over the momentum with the probability distribution ^^{^ (and spin states which is triv- 
ial): 

R^' -n n ^fc''' f ^^Pfc' d^pu> d%, d%i, (0) (0) Rf,i, 

^k'v =9k'gi'-7—r^ / — ^0 ~h' hi' ^k'l' - • [U6 

{27T)^nk' J pl, p\i, p'l, p'\i, Ik'vnk' 

So that to get the mean rate of elastic collisions per particle of the k'-th species with all 
particles in the system one can just integrate the sum of the gain terms in the collision 
integral fl4.22p over ^^^^3^0^^ and average it over spin: 



The Rfn, can be expressed through the integral ( ICISP as 

5e/ (rr\ ,-(0,0,0,0,010,0) 

Rk'i' = 7k'i'Cr{T)nk'ni'Jl,i, 



r3 



9k'9i'7k'i' n — [{zk' - zi,) K2{zk, + zi>) + Zk>zi\zk> + zi,)K3{zk' + zi,)], (D5) 



It represents some sum over all possible collisions. In the case of the same species one factor 7^';/ just 
cancels double counting in momentum states after scattering and another factor 'yk'V also reflects the fact 
that scattering takes place for ("2 ) ~ pairs of undistinguishable particles in a given unit volume. 
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where al^q, is the classical elastic differential constant cross section of scattering of particle 
of the k'-th. species on particles of the /'-th species. For the case of large temperature or 
when both masses are small, Zk' <C 1 and zii <C 1, one has expansion 

Rt'v = 9k'9i'lk'i'^^ [I - ^(4 + 4) + ■■■) ■ (D6) 

For the case of small temperature or when both masses are large, Zk' ^ 1 and zi> ^ 1, one 
has expansion 

~ , V2al\,T^Zk'Z,,/^;7T^e-^>^'-^" ( 84 + l^z^'Zy + 8zf, \ 

^k'v = 9k'gvlk'v i H 7 \ ^ h ... . {iJi) 

TT^I^ V SZk'Zv[Zk' + Zv) J 

For the case when only one mass is large, zii ^ 1, one has somewhat different expansion 

_ ^/^^tvT\zy + l)zTe-^y-^^' ( 44 + 15^,, + 15 \ 
^k'V - 9k'gvlk'v I i H — — z^, ... J . (U8j 



IZk' + 



The a(T)4?/^°'°'°l°-°) in the m 

can be replaced in the limit of high temperatures 

with A'Ka'^iii{\vk'\) and in the limit of low temperatures with 47r(T^',;,(|{4/|) + mk'/rni' = 
^''^'^'k'i'{Wrei,k'i'\), where {\vk'\) is the mean modulus of particle velocity of the k'-th species 



! (PVk' f^\pk') zlK2{zk') y nzk' K2{zk 

and {\vrei,k'i'\) is the mean modulus of the relative velocity, which coincides with the {\vk'\) 
for high temperatures. Then the resultant collision rate Rfj, would reproduce simple nonrel- 
ativistic collision rates know in the kinetic-molecular theory. To get the mean free time one 
has just to invert the Rf, 

4 = 4- (Dio) 
^k' 

The mean free path If, can be obtained after multiplication of it on the 

a = (DID 



For the single- component gas one gets 



jei_ _ ne-^'jzi + l) . 

The nonrelativistic limit of the (ID12p with gi = 1 coincides with the same limit of the 
formula 
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which is the mean free path formula coming from the nonrelativistic kinetic-molecular theory 
obtained by Maxwell. The ultrarelativistic limit of the (]D12p with gi = 1 coincides with the 
same limit of the formula 

If = (D14) 

Analogically one can introduce inelastic rates i?^*^*^' of any inelastic processes to occur for 
the k'-th particles species. Then the mean free time in what any inelastic process for 
particles of the k'-th species occurs, can be introduced as 

= (D15) 

The mean free path for particles of the k'-th species is obtained through the rate Rf, + i?^?^' 
and can be written as 

Del I jDinel ' (D16) 
^k' + ^k' 
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